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ABSTRACT 



We have studied the sensitivity to variations in the the triple alpha and 
12 C(a,7) 16 reaction rates of the production of 26 Al, 44 Ti, and 60 Fe in core- 
collapse supernovae. We used the KEPLER code to model the evolution of 15 M Q , 
20 M©, and 25 M stars to the onset of core collapse and simulated the ensuing 
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supernova explosion using a piston model for the explosion and an explosion en- 
ergy of 1.2 xl0 51 erg. Calculations were performed for the Anders and Grevesse 
(1989) and Lodders (2003) abundances. Over a range of twice the experimental 
uncertainty, a, for each helium-burning rate the production of 26 Al, 60 Fe, and 
their ratio vary by factors of five or more. For some species, similar variations 
were observed for much smaller rate changes, 0.5 a or less. The production of 
44 Ti was less sensitive to changes in the helium-burning rates. Production of all 
three isotopes depended on the solar abundance set used for the initial stellar 
composition. 

Subject headings: gamma-rays: general, nuclear reactions, nucleosynthesis, abun- 
dances, supernovae: general 



Introduction 



Observations of gamma rays from radioactive nuclei, an d, in particular, fr om 26 Al, 44 Ti, 
and 60 Fe, with halflives of 7.2 x 10 5 yr, 60 yr, and 2.62 x 10 6 yr ( IRugel et al.ll2009[ ). respectively, 
provide information about the sites and nature of stellar nucleosynthesis that is difficult to 
obtain in other ways. For example, the spatial distribution of 26 Al inferred from observation 
of its gamma rays indicates that their most probable source is massive stars; a comparison 
with the observed intensity of 60 Fe, 26 Al, and 44 Ti gamma rays provides a stringent test of 
supernova models; a comparison of the observed intensity of the 26 Al gamma rays and the 
calculated production of 26 A1 in core-collapse supernovae provides an estimate of the rate of 
supernovae in the galaxy; the failure to observe 44 Ti gamma rays, other than from Cas A 
may indicate that emission of 44 Ti in supernovae is fairly rare, or that the super n ova rate has 



been ov erestimated. For a discussion of these issues see iLeising fc Diehll ( 120091 ); iDiehl et al. 
(l2006ah . 



Observations have now yielde d accurate gamma ray intensities for these nuclei (IDiehl et al 



2006al and the references therein; IWang et all 120071 ). but the reliability of deductions based 
on these intensities is limited by our understanding of the production mechanisms that 
synthesize 26 Al, 44 Ti, a nd 60 Fe. For examp l e, est i mates of the ratio 6 °Fe/ 2 6 Al have varied 
greatly as discussed by Limongi k, Chieffil (120061 ) ; IWoosley fc Hegerl ( 120071 ) . Early calcu- 



lations by iTimmes et al.l (jl995[ ) are in reasonable agreement with present results, whereas 
some later calculations using the same supernova code yield much different results. A de- 
tailed discussion, t racing these changes to the effects of reaction rate choices, can be found in 
Woosley fc Hegerl ( 120071 ). Additional uncertaint ies arising from dif f erent a strophysical model 
assumptions or approximations are discussed in lWoosley fc Hegerl (120071 ); iLimongi fc Chieffi 
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( 120061 ); lYoung &: Fryerl ( 120071 ). In addition, it appears that y ields may be significant ly dif- 
ferent when simulations are done in two or three dimensions (IMeakin fc Arnettll2006l ). 



The yield uncertainties resulting from reaction rate uncertainties that have been exam- 
ined previously are mostly for reactions directly involved in the nucleosynthesis processes, 
but other reaction rate uncertainti es may also be important . In their recent study of 26 Al 
and 60 Fe yields from massive stars iLimongi fc Chieffil ( 120061 ) showed that the amount of C 
fuel left by the standard helium burning affects all subsequent shell burning episodes and 
the yields of 26 Al and 60 Fe. Since the yields of 12 C and 16 at the completion of core helium 
burning are determined by the rates for the triple-alpha (i^a) and 12 C(a,7) 16 (Ra,i2) re- 
actions, we expect the yields of 26 A1 and 60 Fe to depend on these rates. We have previously 
shown that the production factors fo r many isotopes with A < 90 depend strongly on both 
R 3a and R al2 (TTur et alJl2007l bopgh . 



There have been several studies of the effects of uncertainties in these helium burning 
rates. Only a few have involved changes in R 3a , mostly f or studies of the effec ts of changes 
in coupling constants over long time scales. For example, ISchattl. et al.l (120041 ) investigated 
the effects of changes in the excitation energy of the Hoyle state at 7.65 MeV in 12 C but 
considered only large changes in rate and did not give yields for 60 Fe, 26 Al, and 44 Ti. Other 
authors have dis cussed the results of ch anges in R a} u but yields of 60 Fe, 26 A1, and 44 Ti are 
not often given. Ilmbriani. et al.l (120011 ) evolved a 25 M Q star at two values of R a ,i2 which 
differ by a ratio of 2.35, corresponding to about 0.7 and 1.6 on the rate scales we use for 
R a i2 if one normal izes to the value of t he rate at 300 keV. We will show a comparison with 



the Ilmbriani. et al.l results in Section 3. iHoffman. et al.l (Il999l ) used two different R a ,i2 and 
present results for a 15 M Q star, but also changed many other reaction rates so an accurate 
comparison with our results is not possible. 

It is clear that available studies do not provide a detailed estimate of the variation of 
the yields of 60 Fe, 26 Al, and 44 Ti with changes in the helium burning rates. This shortcoming 
is especially important because one anticipates that the variation of the yields with rates 
might be rapid. The lifetimes of th e long lived isotopes depen d strongly on the temperature 
of their environment. For example, ILimongi fc Chieffil (120061 ). note that the lifetime of 26 Al 
drops to 0.19 yr at log (T/K) = 8.4 and that of 60 Fe to 0.5 yr at log (T/K) = 9. Their yields 
will then depend on whether or not they are formed in a convective region of the star and 
are carried relatively quickly to cooler regions. Since the details of the convective structure 
of a star depend on the helium burning rates, we expect that the production of 60 Fe, 26 A1, 
and 44 Ti might show a similar sensitivity. A detailed study, with closely spaced rate changes 
is necessary. 



In this paper we describe a systematic study of the changes in supernovae synthesis 
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of 26 Al, 44 Ti, and 60 Fe that result from changes in the helium burning reaction rates and 
changes in the initial stellar abundances. The emphasis is not on the absolute yields, but on 
the variations in the yields induced by changes in the helium burning reaction rates and the 
initial abundances. In total about 200 s tars were evolved ; the specific values of the helium 



burning rates used are shown in Fig. 1 of iTur et al.l (120071 ). We find that over the ±2<r range 
of rates, the yield changes are large, often a factor of five or more, and that the yields of 60 Fe, 
and 26 Al often vary rapidly with reaction rate and do not show simple monotonic behavior. 

In §2 we present details of the calculations. In §3, we compare our results for the stan- 
dard values of the helium burning reaction rates to prior calculations (AG89 abundances), 
and then present the outcome of our rate and abundance changes for the three radioactivi- 
ties under study. In §4, we discuss other important uncertainties that can impact the yield 
predictions and compare our results with results of observations. Our conclusions are given 
in 55. 



Stellar models and reaction rates 



We used the KEPLER code rtWeaver et alil978l ; IWoosley fc Weaverlll995l ; iRauscher et al. 



20021 ; IWoosley et al.ll2002l ) to simulate stellar evolution from standard hydrogen burning up 
to core-collapse; a piston placed at the base of the oxygen shell was then used to simulate 
the explosion. The outward velocity of the piston was set to impart a kinetic energy of 
1.2xl0 51 erg to the ejecta as they escape to infinity. After estimating the fallback from 



( 


2007 


). See 


Tur et al. 



(120071 ) for a more detailed description of the models. 



We used the rates of IRauscher fc Thielemann! (120001 ) for the 26 Mg(p, n), 26 Al(n,p), 
26 Al(n, a), 40 Ca(a,7) , 44 Ti(a,7), 44 Ti(a,p), 59 Fe(n,7), and 60 Fe(n, 7) reactio ns, the rates 



of 



Iliadis et all (120011 ) for the 25 Mg(p, 7), and 26 Mg(p, 7) reactions, the rates of iFisker. et al. 
( 20011) for the 44 T i(q,p) reaction, and for 22 Ne(a,n) rate we used the lower limit from 



(j Jaeger et al.ll200ll ). The calculations described in this paper were performed over a signif- 
icant time period. Sometimes the preferred rates changed over that period, but because of 
our desire for a systematic set of studies, focusing on the rates of the helium burning reac- 
tions, we retained our initial choices for the entire set of simulations. We expect that our 
simulations will provide a reasonable estimate of the relative yield uncertainties resulting 
from the uncertainties in the helium burning rates, even if other rates or approximations 



chang e. A good overview over the reaction rates used can also be found in IRauscher et al. 
(120021 ). 
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We calculated 26 Al, 44 Ti and 60 Fe yields for ±2 a va riations of R^a and R a ,u from their 
standard v alues as given by ICaughlan fc Fowler! (119881 ) and by 1.2 times the rate recom- 
mended by iBuchmannl (119961 ). respectively. The ±1 a experimental uncertainties were taken 
to be 12% and 25%. We performed calculations for 15 M , 20 M , a nd 25 stars and 
for three different rate sets: (1) R 3a was kept constant at its value from ICaughlan fc Fowler 



(119881). and Rn 12 wa s varied; (2) R a ,i2 was held constant at 1.2 times the rate recommended 
by IBuchmannl (119961 ) and R 3a was varied; (3) both standard values of the rates were varied 
by the same factor, so their ratio remained constant, to investigate the extent to which only 
the ratio of the two rates is important. 



We also determined the y ields for two different stellar compositions: lLoddersI (120031 ) 
and lAnders Sz Grevessd ( 119891 ) , hereafter L03 a nd A G89. The L03 abundances for C, N, O, 
Ne are rather close to those of lAsplund. et al.l ( 120091 ). 



Results 



3.1. Comparison with previous calculations 



In Table HJ we show the final yields for 2 6 Al, 44 Ti, a nd 60 Fe for our 15 M Q , 20 M , and 



25 M Q stars, and their average using a Scalo ( jScalol 119861 ) Initial Mass Function (IMF) with 
a slope of 7 = —2.6, for the standard values of R 3a and R a ,u, and for both the AG89 and 
the L03 abundances. The difference in the yields due to the initial composition is largest 
for the 25 M star, where the yields are reduced by ~ 28 %, ~ 15%, and ~ 81 % for 26 Al, 
44 Ti, and 60 Fe, respectively, in going from the AG89 to the L03 abundances. This effect is 
significant, but smaller than the effect of variations in the helium burning rates, as we shall 
see in the next subsection. The major difference in the abun dance sets is the substantially 
reduced CNO content of L03 compared to AG89 (see Fig. 2 in lTur et al.ll2007l ). The central 
carbon mass fraction for the 25 M Q star, h owever, changes only by about 5% and is 0.1903 
(0.1790) for the AG89 (L03) abundances ( flur et al.l E)0~7T ) . The IMF-averaged values for 
26 Al and 44 Ti are very similar for the two abundance sets. 

Figure [T] shows a comparison of our results for the yields of 26 Al, 44 Ti, and 60 Fe (for stan- 
dard values of R 3a and R a ,u, and the AG89 abundances), and the result s from some recent 
calcul at ions: in Figure [ p a our 26 Al yields are shown toget h er with those of iLimongi fc Chieffi 



( 120061 ): iRauscher et al. 



J2002I): iThielemann et aH Jl996h : iTimmes et aD dl995h (the latter 
is based on the s urvey by IWooslev &: Weaver! Il995l): in Figure we compare o ur 60 Fe 



Rauscher et al. 



Timmes et al. 



yields to those of Limongi fc Chieffi ( 20061 ): 

in Figure [T]c, the 44 Ti yields are plotted for our study as well as for IRauscher et al 



(J2002|) 



1995h : 
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Thielemann et all (119961 ) . The iThielemann et al.l (119961) yields for Fe are no t shown because 



they w ere much smaller than the yields shown here. iLimongi fe Chieffil ( 120061 ) ; lTimmes et al. 
( 119951 ) only studied the yields for 26 Al and 60 Fe, so no 44 Ti yield can be shown. 



The 26 Al yields found for most of these studies are in rough agreement, given the un- 
certainties due to reaction rate choices and astrophysical model differences described in the 
Introduction. The lar ger differences for 6 Fe c an mainly be explained by the reaction rate 
differences discussed in lWoosley fe Hegerl ( 120071 ) . The yield curves of Figure [TJ for the present 
study have dips at 20 M for the three isotopes studied in this paper. This dip corresponds 
to an inv ersion of the mono t onical ly increasing general trend in the mass-radius relation as 
noted by ILimongi fe Chieffil (120061). A related behavior for the 20 M star was observed for 
medium- weight isotopes ( ITur et al.ll2007l ; iRauscher et al.ll2002[ ). as a result of merging of the 
O, Ne, and C burning shells about a day prior to the supernova explosion. 



3.2. Effect of variations of the helium burning reaction rates 

We ado pt a three-charac ter notation to label our plots, e.g., LA3, similar to the notation 



adopted by ITur et al.l ( 120071 ). The first character can be an "L" (to denote the L03 initial 
abundances) or an "A" (for the AG89 initial abundances). The second character denotes the 
study: "A", when R 3a was kept constant, and R at u was varied, "B", when both rates were 
varied by the same factor, so their ratio remained constant, and "C", when R 0) i2 was held 
constant, and R^ a was varied. The third character is the number of stars included in the 
IMF average; it is typically 3, for averages over the three stars (15 M , 20 M , and 25 M ). 
When no third character is present the values apply to a single star. 

Figures Wfl-f show the yields of 26 Al, 44 Ti, and 60 Fe for the 25 M star as a function of 
reaction rate for the three sets of rates (see Introduction) and for the two abundance sets 
we used: AG89 and L03. Figures |3]a-/ show the same for the IMF- averaged yields. The 
yield often varies non-monotonically and sharply with rate. We find a strong dependence of 
the yields on R 3a and R a ,i2, and a smaller dependence on the initial solar abundances. The 
strongest variations occur for the 60 Fe yields, with weaker variation for the 26 Al yields, and 
still weaker variations for 44 Ti. 

For the L03 abundances, in the LA series for the 25 M star, the 60 Fe yield reaches a 
maximum value of 3.59OxlO _4 M when R a ^2 is reduced by 24% (relative to its standard 
value), and a minimum of 8.285 xlO _6 M when it is reduced by 48%, a factor of ~ 43 
change in yield. Changes in R 3a can cause similarly large effects: in the LC series for the 
25 M star, the 60 Fe yield is at it maximum value of 2.400x10 4 M when R 3a is increased 
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by 24% (relative to its standard value), and a minimum value of 2.881 x 10 -5 M when R 3a 
is reduced by 6 %, a factor of ~ 12 change. Even for 44 Ti, differences of a factor of about 
2 are observed for certain cases. The differences for the IMF averaged case are somewhat 
smaller but still substantial, especially for 60 Fe. 

For the AB and LB series, where both rates are a given multiple of their standard 
values, the variations are weaker. The variation of the central carbon abundances with 
reacti on rate is about one-third of that for comparable variation of the reaction rate R a ,i2 



alone (ITur et al.l (120071 )). This is consistent with the generally weaker variations in the yields 



of the radioactive nuclides. 

Since the half-lives of 26 Al and 60 Fe (~ 10 6 yr) are much longer than the average time 
between two galactic SNII events (~ 2 per century), a steady state assumption, where the 
production rate balances the decay rate is reasonable. The changes in the production yields 
discussed above then translate directly into changes in the predicted photon flux. 

More detailed results for variation of R% a and R a ,i2, and for the two initial abundance 
sets are given in Tables fS] and 



We can compare our results to those of llmbriani. et al.l (120011 ) at R at u = 0.7 and 1.6 



Since we did not do simulations at these exact values, we have read values from Fig. [2]a. 
The results, especially for 60 Fe, are subject to large extrapolation uncertainties. The energy 
dependencies of the reaction rates also differ. The ratio of our results at 0.7 to those at 1.6 



for 26 A1 ( 60 Fe) is 0.8 (0.6) and those of llmbriani. et al.l is 0.17 (1.22). Individual values of 



the yields agree within a factor of 2 — 3. It is difficult to judge whether these differences are 
significant, given the different assumptions on stellar physics that entered the calculations, 
as outlined in the introduction. 



3.3. Production mechanisms 

Figure H] shows the yields of 26 Al, 44 Ti, and 60 Fe for a 25 M Q star using the L03 abun- 
dances at various stages of stellar evolution for the standard values of Rs a and R a ,u (Fig- 
ure©) and when R 3a is increased by 18 %, R a ,i2 maintained at its standard value (Figure©). 
The stellar burning stages for which yields are plotted and that are also used in Figure \5\ 
are: 

He-ign: At this time 1 % of the helium has burnt in the center. 

He-dep: Just before central He depletion (when the central He mass fraction has reached 1 %). 
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C-ign: Just before central C ignition (when the central temperature reaches a value of 5 x 10 8 K) . 
Central helium burning has finished. 

C-dep: Central C depletion (when the central temperature has reached 1.2xl0 9 K). Central 
C burning contributes mostly to those regions of the star below the mass cut. 

O-dep: At central O depletion (when the central O mass fraction has dropped below 5%). 

Si-dep: Central Si depletion (when the central Si mass fraction drops below 10 -4 ). 

pre-SN: At the pre-supernova stage when the contraction speed in the iron core reaches 1000 km s^ 1 . 
This is typically about 0.25 s before core bounce. 

57V unmix: 100 s after the passage of the shock wave. At this time all thermonuclear reactions 
and neutrino-process have essentially ceased, but Rayleigh- Taylor (RT) instabilities 
that cause mixing after the passage of the shock have not yet mixed the ejecta (in our 
approximation). This is only shown in Figures [5]g?, [5]e. 

Final: > 100 s after the pas sage of the shock way e. This includes mixing due to RT instabilitie s 
parametrized as in (IRauscher et al.ll2002l ) and fallback (see also lYoung &: Fryerl 120071 ). 
In Figures \5)p,f this has been labeled as "S7V mix" . 



The convective history for the above two cases is shown in Figure \5)p,, b and the times of 
the snap-shots are indicated. Figure 0c, d shows the pre-supernova structure for these two 
cases including convective regions and mass fractions of 26 Al, 44 Ti, and 60 Fe. The inner dark 
gray region indicates iron core where we no longer follow detailed nucleosynthesis as these 
regions become part of the collapsing iron core and disappear in the remnant. Figure E]e,/ 
shows the distribution of the isotopes inside the star at the different evolution stages. Here, 
light gray regions indicate where detailed nucleosynthesis was no longer followed. Notably, 
in both of the last two rows of Figure |5] the vertical solid black line indicates the location of 
the "mass cut" after the supernova explosion - everything below that line disappears in the 
remnant. The solid lines and symbols in Figure H] show only the contributions outside the 
this mass cut, the dashed line and hollow symbols show mass contained in the entire star as 
computed at each stage. In general, the piston in the SN explosion can be below the mass 
cut if there is "fallback" - mass that has first moved outward during the SN explosion, but 
then did not escape and fell back onto the compact remnant instead; in the specific cases 
shown here there was not such "fallback" . 

As is shown in Figure |5]a, b changes in reaction rates can cause significant differences in 
the convective history of a star. Convective zones can carry radioactive nuclei to cooler zones 
of the star where their effective lifetime is longer and their survival is enhanced. It can bring 
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nuclei to hot regions where they can be destroyed, or transport "fuel" to very hot regions to 
synthesize nuclei. In our case we find a convective oxygen-burning shell visible at roughly 
2 M . In Figure 0a, this shell starts at approximately log(time till core collapse / yrj^J of 
—3 and mixes with the carbon-rich layer above, whereas in Figure Eh where R% a is increased 
by 18 % this burning layer starts at log(time till core collapse / yr) of —3.5 but does not mix 
with the carbon-rich layer above. In the case of merged layers 60 Fe is efficiently destroyed 
whereas it remains unchanged in the separated carbon layer (Figures |5]c and EH). 44 Ti, on 
the other hand, is efficiently produced at the bottom of the hot merged layer and mixed out 
throughout the entire convective domain, whereas in the case of separated layers only some 
of it is made at the bottom of the oxygen-burning layer and little is made in the carbon 
layer. The merging of the layers is also not beneficial to the production of 26 Al; this works 
better at the bottom of the separated carbon-burning layer of the star with the increased 

Rsa- 



3.4. 



2(> 



Al 



A substantial amount of 26 Al is already present at core helium ignition as a result of 
the 25 Mg(p, 7) 26 A1 reaction occurring during core H burning. This reaction is essentially the 
only means of production of 26 Al, even during the later stages of stellar evolution. During 
core helium burning we see most of 26 A1 in the H envelope with an abundance peak at the 
location of the hydrogen-burning shell. In the helium-burning core 26 Al is mostly destroyed. 
The amount of 26 Al remains essentially constant until core O depletion and then increases 
in the later stages of stellar evolution where 26 A1 is made in C shell burning. The amount 
finally ejected after explosion is about an order of magnitude larger than the amount present 
at helium ignitio n. A discussion of the reac tions that create and destroy 26 Al in later burning 
stages is given in lLimongi fc Chieffil (120061 ) . 



3.5. 60 Fe 

Only a negligible amount of 60 Fe (~ 3xlO~ 15 M ) is present at core helium ignition. 
That amount rapidly increases to ~ 1 x 10~ 7 M at core helium depletion and gradually builds 
up through the later stages until explosion. Most of it is found at the lower part of the split 
He-burning shell and is made after central oxygen depletion. This last point underscores the 



^^Here we assume that the core collapses in 250 ms after the end of our pre-SN calculation, when the iron 
core reaches an infall velocity of 1000 km/sec 
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importance of virtually every stage in the star's life for the production of this isotope. 60 Fe is 
mainly produced by neutron capture on 59 Fe and destroyed by 60 Ni(n, 7)reactions. When the 
standard values of the helium burning rates are used, the pre-SN results show the destruction 
of some 60 Fe; when R 3a is increased by 18 %, this does not occur. This is due to the merging 
of the oxygen-burning and carbon-burning layers in the standard case which is not present 
in the case with the increased R^ a (Figure [5]). 60 Fe from the entire carbon-rich convective 
zone can be mixed to the bottom of the very hot oxygen layer where it is partially destroyed. 
The 60 Fe present in the He shell is not directly affected by this destruction, but due to the 
merging layers, the shell gets less hot and hence produces less 60 Fe (Figure Oe,/). There is a 
significant difference between the amount of 60 Fe ejected in those two cases: 3.933 xlO _5 M 
and 2.201 xlO -4 M , respectively, just as a result of variations in the helium burning rates. 



3.6. 44 Ti 

The amount of 44 Ti present at core helium ignition is minuscule - ~ 0.25 g. It increases 
through the later burning stages, but remains below lxlO _8 M Q outside of the mass cut 
until core Si depletion for the case R^ a = 1.18, and still later for the standard values of the 
helium burning rates. 44 Ti is only made by O burning and later phases that all start off 
below the SN mass cut . A dis cussion of the reactions that create and destroy 44 Ti is given in 



Vockenhuber. C. et al.l ( 120081 ). Remarkably, there is a significant amount of 44 Ti at the pre- 
SN stage just outside of the Fe-core in the Si-Shell. All of this disappears inside the remnant 
in these two cases; only some 44 Ti that is made in O shell burning is being ejected. This 
contribution is significantly bigger in the case of the merged O/Ne/C layers of the standard 
case as compared to the split layers of the case with R Sa = 1.18. Whereas the standard 
case has lOx more 44 Ti in the pre-SN model outside the SN mass cut than the R% a = 1.18 
case, the final production in both cases is entirely dominated by SN explosions. It would, 
therefore, be sensitive to uncertainties in the explosion energy and the asymmetries in the 
SN, and how these depend on a pre-SN structure, which is quite different in the two cases. 
That is, if the explosion was different in the two cases and not a "standard" 1.2xl0 51 erg 
explosion as we have assumed here, the 44 Ti yields might be quite different. 



4. Discussion 

We have shown that the uncertainties in the helium burning reaction rates cause large 
differences in the predicted yields of 26 Al, 44 Ti, and 60 Fe, with the changes in 60 Fe being 
particularly large. These yields are also sensitive to uncertainties in the rates of other im- 
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portant reactions and i n the approximate treatment of convection (ILimongi fc Chieffil 12006 
Woosley fc Hegerl 120071 ); there may also be sites other than massive stars that contribute to 
the abundances of 26 Al, 44 Ti, and 60 Fe. The changes we have examined here, however, give 
a lower limit on the overall uncertainties which is sufficiently large so as to limit deductions 
that depend on their yields. Including other uncertainties will only make the situation worse. 

We considered only stars with masses less than 3 M m and consequently di d not include 
the contributions of Wolf-Rayet stars. In the models of lLimongi fc Chieffil (12006). Wolf-Rayet 
winds contribute less than half the yield of 26 Al in these heavier stars; IWoosley &: Heger 
( 120071 ) find a similar contribution of wind and explosion for a 60 M Q star. It seems reasonable 
to assume that the explosive process in Wolf-Rayet stars is subject to uncertainties similar 
to those we found for lighter stars, and that their overall contribution will be subject to 
similar uncertainties. This remains to be confirmed by further calculations. 

Finally we consider two comparisons of the nucleosynthesis yields and astrophysical 
observations. As was mentioned above, in the steady state expected for isotope lifetimes short 
compared to galactic evolution times, the production r ate of gamma rays is proportional to 
the amount of material made in a supernova explosion. iDiehl et al.l ( 12006bl ) have shown that 
the observed 26 A1 gamma ray flux corresponds to a supernova rate of 1.9 ± 1.1 events/century 
for an assumed value of the supernova yield. Since the gamma ray flux is proportional to 
the supernova yield for 26 Al, it is subject to the uncertainties in the yield shown in Figure 
Wfl-f'i we must conclude that the uncertainty in this estimate of the supernova rate is much 
too small. 

A second comparison is to the ratio of the gamma fluxes from 26 Al and 60 Fe. If one 
assumes that for both nuclei the gamma rays come from supernovae, and that the contribut- 
ing supernovae have the same distribution in space, then the ratio of the fluxes is the ratio 
of the numbers of nuclei produced, or flux(60)/flux(26) = (26/60)(yield(60)/yield(26)). The 
yield ratios are shown in Figure [6]a- d. These rati o s sho uld be compared to 60/26 times the 
flux ratios summarized in Table l2lof I Wang et al.l ( 120071 ). namely (60/26)(0.15 ± 0.06). The 
ratios shown in Fig. [6] are much larger for most values of the helium burning reaction rates, 
although the smallest values of the 12 C(a,7) 16 reaction rates do give ratios approaching 
the observed values. But given the size of the variations shown, it will not be possible to 
show that this value agrees with supernovae predictions until the helium burning reaction 
rates are much better known. 



A different sort of comparison arises for 44 Ti. It is seen in the supernova remnant 
Cas A, bu t not elsewhere. Given conventional supernova rates, additional occurrences were 
expected (IThe et al.ll2006l ). 44 Ti production is not particularly sensitive to the helium burn- 
ing rates-see Figure [3]a-/. Possible explanations include a lower production rate as occurs for 
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some reaction rate combinations, or a lower supernova rate than assumed. As was pointed 
out above, 44 Ti produ ction may be particularly sensitive to the supernova explosion physics. 



Young fc Fryer! (120071 ) have shown that these effects can be very large. A reference to our 



explosion model including supernova burning te mperatures and time - scales as well as sensi- 



tivity of yields to explosion energy can be found in lWoosley &: Weaver! ( 119951 1 ; lRauscher et al. 
(l2002h . 



5. Conclusions 

We explored changes in the core-collapse supernovae yields of 26 Al, 44 Ti, and 60 Fe arising 
from variations in the triple alpha and 12 C(a,7) 16 reaction rates. We find that these 
changes are significant, sometimes more than an order of magnitude from the minimum to the 
maximum yield within the ±2 a range of rate uncertainty considered here for both reactions. 
The changes were most important for 60 Fe, smaller for 26 Al, and still smaller for 44 Ti. The 
yield-rate relationship was found to be complex, i.e., non-monotonic, and sometimes with 
large changes for small changes in reaction rates. For this reason it is difficult to give a 
simple statistical uncertainty for the effects of changes in a reaction rate. The high and low 
predictions provide a reasonable estimate for a given uncertainty range. For the triple alpha 
reaction la rate changes are 1.0 ± 0.12 and for the 12 C(a,7) 16 reaction, 1.2 ± 0.3. Thus 
for a 25 M star, the 60 Fe yield varies by a factor of 5.0 (3.4) for ±1 o variations of the triple 
alpha ( 12 C(a, 7) 16 0) rate. Other ranges can be determined from Tables [2] and [3j Obviously 
the detailed dependencies of the isotope production on the rates we varied may differ in 
other model calculations as they also depend on the stellar physics and modeling details. 

We have isolated the effects on 26 A1, 44 Ti, and 60 Fe yields, of variations in the helium 
burning reaction rates within their experimental uncertainties. These effects turn out to be 
large and limit the conclusions we can draw from comparisons with observational data. They 
provide a lower limit to the total uncertainties that will remain even if all the uncertainties 
mentioned earlier are eliminated. This highlights the pressing need for improvements in our 
knowledge of the triple alpha and 12 C(a,7) 16 reaction rates. 

We also explored the dependence of the 26 Al, 44 Ti, and 60 Fe yields on the initi al solar 



abundances us ed for the initial composition of the star: lAnders &: Grevessd ( 119891 1 versus 



Loddersl (120031 ) . There again the dependence is complex, but smaller than that induced by 



changes in the helium burning reaction rates. 



We thank Robert Hoffman for providing the solar abundance files used in this study 
and Stan Woosley for helpful discussions, including studies on the relative influence of the 
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Table 1. Yields (M )xlO 5 : standard values of R^ a and R a ,i2- 



15 Mr., 20 Mp) 25 Mm 3 star ave. 



AG89 abundances: 



26 Al 


2.789 


2.934 


9.781 


4.349 


44 Ti 


4.221 


3.333 


5.707 


4.266 


60p e 


6.333 


4.748 


20.76 


8.965 






L03 abundances: 




26 Al 


2.443 


5.962 


7.023 


4.535 


44 Ti 


4.263 


4.039 


4.846 


4.322 


60p e 


7.481 


1.727 


3.933 


4.924 



a IMF-averaged over the 15, 20, and 25 M Q models. 



Table 2. Yields (M ) for 26 Al, 44 Ti and 60 Fe as a function of the R 3a and R a ,i 2 rates for stars of 15 M , 20 M , and 

25 M n with AG89 initial abundances 



15 Mq 20 Mq 25 M Three Star Average 



Rate Multiplier 51 
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■05 
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■05 
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3.79E 


■05 
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■05 


3 
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■05 
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■01 


6.29E+00 



a The rate multipliers are given in the format (R3 a multiplier, Ti a ,i2 multiplier), where the rates being multiplied are the standard rate values chosen as described in the text 



Table 3. Yields (M ) for 26 Al, 44 Ti and 60 Fe as a function of the R 3 „ and R a ,i 2 rates for stars of 15 M , 20 M , and 

25 M with L03 initial abundances 



15 Mq 20 Mq 25 M Three Star Average 
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a The rate multipliers are given in the format (R3 a multiplier, Ti a ,i2 multiplier), where the rates being multiplied are the standard rate values chosen as described in the text 



- 18 - 




2.5 




■■■■><■" 
--Q- 


This paper 

Rouscher et ol. (2002) 
Thielemann et al. (1996) 


44 Ti 


(c) : 


2.0 












1.5 












1.0 












0.5 




H*' 








0.0 




x 






*5 ; 



14 16 18 20 22 24 26 

stellar moss (M ) 



Fig. 1. — Yields versus stellar mass: comparison between the present study (standard values 
of the helium burning reaction rates, AG89 abundances) and previous publications, (a) 26 Al. 
(b) 60 Fe. (c) 44 Ti. 
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Fig. 2. — Yields versus reaction rate for the 25 M model. For details see the text, (a) AA 
series, (b) LA series, (c) AC series, (d) LC series, (e) AB series, (f) LB series. 
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Fig. 3. — Yields versus reaction rate, averaged over 3 stars (15 M Q , 20 M Q , and 25 M Q ). 
details see the text, (a) AA series, (b) LA series, (c) AC series, (d) LC series, (e) 
series, (f) LB series. 



For 
AB 




He-ign He-dep C-ign C-dep O-dep Sl-dep pre-SN Final 
stages 




He-lgn He-dep C-ign C-dep O-dep Si-dep pre-SN Final 
stages 



Fig. 4. — Solid symbols and lines give yields (total mass of the isotope outside the final mass 
cut of the SN, see Figure |5]) at various evolutionary stages for a 25 M star using the L03 
initial abundances. The definition of the stages labeled on the abscissas are given in the 
text. The hollow symbols and lines show the total yield inside the star (including wind) as 
shown in the bottom row of Figure (a) Standard values of the helium burning reaction 
rates (i.e., multipliers for R 3a and R a ,n'- 1-0 and 1.2, respectively), (b) R 3a higher by 18% 
(i.e., multipliers for R 3a and R a ^2'- 1.18 and 1.2, respectively). 
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Fig. 5. — Evolution of and isotope distribution in 25 M stars of L03 abundances. Left 
Column: Standard values of the helium burning reaction rates, (R3 a ,R a ,i2) = (1.0,1.2). 
Right Column: R 3a enhanced by 18%: (R3 a ,R a ,i2) = (1.18,1.2). Top Row: Convective 
history as a function of time until core collapse. The ordinate is the mass coordinate from the 
center of the star. The green hatched fully convective, and the red cross hatched 

semiconvective. The blue and purple shading indicate net energy generation from 
burning and n e utrino losses, with blue positive and purple negative. For more details see 



Woosley et al.l ( 120021 ). The vertical lines show the snapshots we use and are explained in 



more detail in the test. Middle Row: Stellar structure at the pre-supernova stage as a 
function of the enclosed mass (truncated at 10 M Q ). Dark gray indicates the neutronized 
regions in which nucleosynthesis was no longer followed. The vertical black line shows the 
final mass cut of the SN. The light gray indicates convective regions. The solid gray line 
shows the nuclear energy generation and the dashed gray line shows the density. Bottom 
Row: Distribution of isotopes inside the star as a function of mass coordinate for different 
evolution stages (y-axis, time increased downward). Gray shading indicates neutronized 
regions (until pre-SN) or location of the SN piston (SN lines) in which nucleosynthesis was 
no longer followed as those regions become part of the remnant. The solid black line shows 
the final mass cut of the SN. For the two stars shown here there was no fallback, so the mass 
cut and the piston mass coincide; this is not the case in general. 
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Fig. 6. — Mass ratio of yields of 60 Fe and 26 Al versus reaction rate, averaged over three stars 
(15 M , 20 M , and 25 M Q ). For details see the text, (a) AA series, (b) LA series, (c) AC 
series, (d) LC series. 



